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Abstract 

The  migration  of  slab,  cylindrical  and  spherical  zones 

through  a  block  of  solid  under  the  influence  of  a  temperature  gradient 

has  been  analyzed  theoretically.  The  migration  rate  is  found  to  be 

directly  proportional  to  the  temperature  gradient  in  the  system  and  to 

* 

depend  upon  orientation  -  dependent  atomic  kinetics  of  the  melting  and 
freezing  processes.  She  magnitude  of  the  atomic  kinetic  effect  depends 
upon  the  size  of  the  molten  zone.  Stable  zone  migration  occurs  only  in 
certain  crystallographic  orientations.  The  theory  is  discussed  in  terms 
of  the  limited  experimental  results  available. 


THE  MIGRATION  OF  A  LIQUID  ZONE  THROUGH  A  SOLID:  Part  I 

■  W.  A.  Tiller 

Since  the  initial  paper  "by  Pfann^  on  "Temperature  Gradient 

Zone  Melting"  (TGZM),  only  a  few  publications  have  appeared  utilizing  or 

(2)(3)(4)  (2)(3) 

studying  this  technique..  However,  the  work  of  Wernick'  '  '  has 

shorn  that  the  migration  rate  of  the  zone  is  not  primarily  diffusion  controlled 

hut  depends  upon  the  zone-size .and  zone-orientation  as  well.  Further,  the 

(5) 

work  of  Mlavsky  and  Weinstein'  '  and  others  indicate  that  this  might  he  a 
useful  technique  for  preparing  small  crystals  of  certain  complex  materials. 

The  author  feels  that  the  migration  of  a  liquid  zone  through  a 
solid  under  either  a  thermal  field  or  an  electric  field  driving  force  can 
he  developed  into  a  very  potent  technique  for  both  the  preparation  of 
crystals  and  the  determination  of  fundamental  information.  Therefore,  the 
present  paper  is  an  attempt  to  evaluate  the  scope  of  the  technique  through 
a  theoretical  analysis  of  zone  migration  in  a  temperature  gradient;  in  a 
companion  paper, ^  the  electric  field  effects  are  treated. 

Theory 

The  starting  geometrical  and  thermal  configurations  for  a  slab 
zone  in  a  block  of  solid  may  he  as  illustrated  in  Figure  1.  The  upper 
portion  of  the  block  consists  of  solid  solvent  and  perhaps  some  solute. 

The  zone  region  consists  of  a  substantial  alloy  of  solute  and  solvent  and 
the  lower  block  usually  consists  of  a  single  crystal  seed  of  the  solvent. 


J 

I 


The  ensemble  is  placed  in  a  temperature  gradient  such  that  the  temperature 
of  the  zone  region  is  above  the  lowest  melting  temperature,  T  ,  of  the 
system,  and  such  that  the  highest  temperature,  T.  ,  in  the  system  is  below 
the  melting  point  of  the  upper  or  lower  blocks.  Under  these  conditions, 
the  zone  region  melts  and  the  zone  adjusts  its  length  until  the  zone  con¬ 
centration  is  between  and  (Figure  l).  Reactions  will  occur  at  the 
two  solid-liquid  interfaces  to  produce  the  equilibrium  concentrations 
and  Cg  at  and  T^  respectively.  A  solvent  concentration  gradient  is  set 
up  in  the  zone  and  the  combination  of  solution  reaction  at  T»,  diffusion 
in  the  zone  and  freezing  at  T^  causes  the  molten  layer  to  migrate  through 
the  block. 

In  this  treatment,  the  direction  of  the  maximum  temperature 
gradient  will  be  taken  as  the  Z  axis.  The  migration  of  slab,  rod,  and 
spherical  zones,  having  dimensions  J?.  cms,  as  illustrated  in  Figure  2  will 
be  considered.  The  zones  are  considered  to  have  a  concentration  between 
and  and  to  be  located  at  Z  =  Zq  at  time  t  =  tQ.  The  assumptions  to 
be  made  in  the  treatment  are:  (i)  we  are  dealing  with  a  binary  or  pseudo¬ 
binary  system,  (ii)  no  exchange  of  either  solute  or  solvent  exists  between 
any  portion  of  the  sandwich  and  the  gas  phase,  and  (iii)  the  zone  migration 
occurs  at  such  a  slow  rate  that  the  influence  of  the  interface  motion  on 
diffusion  in  the  zone  can  be  neglected.  Let  us  first  consider  the  case 
where  the  atomic  kinetics  are  isotropic  so  that  migration  occurs  in  the 
0=o  direction. 


1.  Isotropic  Atomic  Kinetics 


The  motion  of  the  zone  occurs  as  a  result  of  (a)  interface 
reaction  to  produce  local  equilibrium  and  (b)  solvent  transport  across 
the  zone  as  a  result  of  the  solvent  distribution  produced  by  (a).  Thus,, 
expressions  need  to  be  developed  only  for  the  rate  of  solvent  transport 
across  the  zone,  this  rate  being  properly  correlated  with  the  rates  of 
melting  and  freezing  at  the  interfaces.  Further,  the  atomic  kinetics  of 
these  interface  reactions  relating  the  reaction  rate  to  the  departure  from 
equilibrium  must  be  considered.  We  shall  treat  the  three  cases"  (i)  uniform 
interface  motion,  (ii)  layer  formation  by  screw  dislocations,  and  (iii) 
layer  formation  by  two-dimensional  nucleation.  For  these  mechanisms,  the 
atomic  kinetics  are  described  by  the  following  equations  (7) (8) (9) 

(i)  V  =  p(T*  -  T)  =  m  p(C*  -  C)  (l) 

where  V  is  the  freezing  velocity  in  the  Z  direction,  T*  and  T  are  the 
equilibrium  and  actual  temperatures  at  the  zone  face,  p  is  the  atomic 
kinetic  coefficient  (l  <  p  <10^  cm/sec/°C  usually),  m  is  the  liquidus 
slope  for  the  solute  (m  is  negative)  and  C  and  C*  are  the  equilibrium  and 
actual  solvent  concentrations  respectively  at  the  zone  face. 

(ii)  V  =  p(T*  -  T)2  =  m2  p(C*  -  C)2  (2) 

-5  -1 

with  p  negative  for  melting  and  with  10  J  <  j  p  j  <  10  generally  and  for 
a  diffuse  layer  edge,  |p|  maybe  as  large  as  unity. 


(iii)  V  =  u  e 


(3) 


-A/(T*  -  T).  _  -A/m  (C*  -  C) 


with  p  and  A  negative  for  melting  and  where  |  p  |  is  the  product  of  the 
interfacial  area,  the  layer  step  height  and  the  pre- exponential  rate 
term  for  two-dimensional  nucleation.  (10  <  j  p  |  <  10 ^  usually)  and  [  A  j 
is  given  by 


A 


S  A  <T  2  V  T* 
u  e  m 

k  T  A  H 


(4) 


In  eq.  k,  S  is  a  geometrical  factor  (S  ~  Tt),  A  is  the  layer  step  height, 
Q-  is  the  layer  edge  surface  energy,  is  the  specific  volume,  k  is 
Boltzmann's  constant  and  AH  is  the  latent  heat  of  fusion  (l  <  J  A (  <  10 ^ 
usually).  In  Figure  3,  the  general  equilibrium  temperature-concentration 
relationship  is  given  for  the  case  where  atomic  kinetics  have  been  con¬ 
sidered. 

Since  the  migration  of  the  zone  occurs  at  such  a  slow  rate,  both 
the  temperature  and  solute  fields  in  the  zone  are  very  well  described,  at 
any  moment  of  time,  by  the  solution  to  Laplace's  equation.  Using  straight¬ 
forward  procedures,  the  temperature  for  all  three  zone  shapes  is  found  to 
depend  only  on  Z.  Assuming  the  same  thermal  conductivity  outside  the  zone 
and  constant  parameters  in  the  two  phases,  the  temperature  gradient  in 
the  zone,  G,  may  be  related  to  the  temperature  gradient  in  the  bulk  solid, 
G*. 


K  V  AH 

G  -  jr  G*  -  -4- — 

KI.  ■ 


(i) 


Slab 


(5) 
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(ii) 


Cylinder 


/ 

(6)  . 


3K„  G*  VK  A  H 

(iii)  Sphere  G  =  ^  (7) 

Where  is  the  velocity  of  the  freezing  face,  is  the  thermal  conduc¬ 

tivity  of  the  zone  and  Ka  is  the  thermal  conductivity  of  the  solid.  From 

-  O 

Figure  k,  ve  can  see  that,  in  the  velocity  range  of  our  interest,  G  is 
practically  unaffected  by  the  magnitude  of  (V^  <  10  ^  cm/sec). 

The  solvent  concentration,  C,  may  be  readily  shown  to  depend 
only  on  Z  (since  the  boundary  conditions  depend  only  on  Z)  for  all  zone 
geometries.  This  is 

C  =  A  Z  +  B  (8) 

i 

where  A  and  B  are  constants  that  vary  slowly  with  time.  Let  C  represent 
the  initial  solvent  concentration  throughout  the  upper  block  of  solid  and 
let  subscripts  refer  to  the  value  of  0  under  consideration  where  ©  refers 
to  the  outward  normal  to  the  zone,  i.e.  for  the  slab  zone,  0  =  0  and 
0  =  7t?  refer  to  the  upper  and  lower  faces  respectively.  The  conservation 
of  solute  conditions  give 

V-,X  d  -  V  °« -  -D  <!>*  -  7  K  -  C-K>  <9) 

<°o  -  cl>  -  J  <°o  -  °K> 


(10) 
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where  kQ  is  the  phase  diagram  solvent  partition  coefficient  (kQ  >  l)  and 

D  is  the  solvent  diffusion  coefficient  in  the  melt.  We  can  see  immediately 

from  eqs.  9  and  10  that,  since  Cq  and  are  functions  of  time,  VQ  and 

will  also  he  time  dependent.  Further,  since  VQ  and  will  generally 

he  slightly  different  in  magnitude,  X  will  he  time  dependent  also.  The 

ratio  of  V_,  /V  is 
rC  '  o 

Z2L .  -i  (c°  -  cl)  ' 

vo  “trrsj  c„  (11) 

In  Figure  5,  the  ratio  -V^/Vq  is  plotted  versus  (C*  -  Cq ) /C ^  for  several 


values  of  k  . 

o 


From  eqs.  1,  2,  and  3  written  in  terms  of  and  Vq,  both  Cq 
and  can  he  determined,  i.e.  Vq  =  m  |iQ  (Cq  -  CQ),  V^.  =  m  p.  ^  (C-^  -  C^), 
Inserting  Cq  -  into  eq.  9,  we  find  the  following  relationships 

(i)  Uniform  interface  advance 


L£  G/m  ] 


1  ‘  "  k[ 


(ii)  Screw  dislocation 


V  * 


l1  -  ^ 


-  [  €  G/m] 

(v,r,r1/2  <- 


♦  iAS) 


(iii)  Two-dimensional  nucleation 


1  -  k  < V"1 


-■  [€G/ffi] 
A 


TV  TC; 


I 

I 

I 

[ 


0 


where  0  -  (TQ  -  T*)//  ,  €  -  -D/(l  -  kQ)  C  ^  and  ^  -  \iQ  VW/VQ.  Here, 

€.  is  almost  a  constant  since  C  7^  «  C*_,  which  is  not  likely  to  change  hy 
more  than  ten  per  cent  in  a  single  experiment.  In  general,  melting 
kinetics  are  much  faster  than  freezing  kinetics  and  we  can  often  neglect 
terms  containing  relative  to  terms  containing  p  ^ .  Since  C  ^  varies 
slightly  over  the  curved  zones,  these  zones  will  become  slightly  elliptical, 
with  the  Z-axis  as  the  semiminor  axis  of  the  ellipse.  From  eqs.  12,  13, 
and  14,  we  can  see  that  the  effect  of  the  atomic  kinetics  is  to  slow  down 
the  zone  motion  and  make  a  function  of  zone  size,  J> .  In  Figure  6, 
eqs.  12  and  13  have  been  plotted  in  the. form  of  versus  J;  for  several 
values  of  the  atomic  kinetic  parameter  P  =  (l/p7r  -  l/p'  )  for  eq.  12  and 
P  =  (1/^/iTp  +  1/^ySJ  )  for  eq.  13.  We  see  that  for  uniform  interface 
kinetics,  V7f.  is  linear  with  X  for  small  X ,  whereas  for  screw-dislocation 
kinetics,  is  sigmoidal  with  Jb  .  In  Figure  7,  eq.  14  has  been  plotted 
as  versus  JL  for  several  values  of  and  p  T(.  assuming  that 
-Ao/ln  (p’/VTC)«  A^/ln  (p^/V-^).  Here  we  see  an  even  greater  si©noidal 
relationship  between  V^-  and-X  and  that  the  atomic  kinetic  effects  are 
observed  at  an  order  of  ' magnitude  larger  values  of  X  . 

In  general,  the  zone  length  will  vary  with  time  and,  utilizing 
ah  overall  conservation  of  solvent  condition,  we  can  develop  an  approxi¬ 
mate  expression  for  .  Denoting  X{o)  and  JL{t)  as  the  zone  lengths  at 

ft  (v o  +  V  )  dt 

t  =  o  and  at  t  and  the  zone  travel  distance  Z(t)  =  J  q - — - 

( z  (t)  [c1  -  kQ  (cM  (0)  -  G*  Z  (t)/an)J  +i(o)  cM  (0) 

|  4(0)  G*  Z  (t)/i 


Mt)  - 


m 


,  we  find 

[(15) 


) 
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where  (0)  is  the  initial  concentration  of  the  mid-point  of  the  zone. 
Since  /  (t)  refers  to  the  dimensions  in  the  Z-direction,  the  circular 
cylinder  will  become  elliptical  and  the  sphere  ellipsoidal  due  to  changes 
in  ,.(<■  with  time. 


2.  Convection  Effects 


The  foregoing  treatment  has  dealt  with  the  case  of  zero  con¬ 
vective  mixing  in  the  zone.  If  the  convective  mixing  were  complete,  then 
the  concentration  would  be  constant  throughout  the  zone.  In  this  case, 
we  have  for  uniform  interface  advance  kinetics 


and 


G  JL 


_v  _  v  «  MM. 
o  at. 


(16) 


(17) 


For  the  case  of  partial  mixing,  we  must  replaced  by  a  &c  in  eqs.  9  and 
10  where  a  ?=^l/2  and  Sc  is  the  convectionless  layer  thickness  defined  in 
Figure  8.  In  this  case,  eq.  12  becomes 


-  [  €G/m]  2  X/Sc 

M  =  r  2  f.  ,  1  ITT 


(JL.  .  A.) 

V,,  .  M  *  / 


(18) 


m6c  K 

t  From  eq.  18,  we  can  see  that  convective  mixing  enhances  and  may  be 
extremely  helpful  in  the  determination  of  the  atomic  kinetic  coefficients 
for  the  cases  there  -€/im  (l^cc  1  by  making  5c/l£«  1.  For  the  other 
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atomic  mechanisms,  the  appropriate  equations  replacing  eqs.  13  and  l4  may 
he  readily  derived  by  substituting  Sc/2  for  JL  in  the  denominator  and 
multiplying  by  2 JL/&  in  the  numerator. 

3»  Anisotropic  Atomic  Kinetics 

Our  interest  in  the  effects  of  anisotropic  atomic  kinetics  may 
be  directed'  in  two  directions,  (i)  do  they  introduce  microscopic  shape 
changes  into  the  zone  contour  and  (ii)  do  they  cause  macroscopic  shape 
changes  in  the. zone  contour. 

When  the  atomic  kinetics  are  isotropic,  microscopic  undulations 
in  the  zone  surface  contour  are  unlikely  to  occur  since  any  changes  in 
curvature,  of  either  sign,  will  alter  the  equilibrium  temperature  in  such 
a  way  as  to  force  the  surface  back  to  its  original  contour.  Further,  no 
macroscopic  shape  changes  can  occur  since  the  fastest  melting  direction 
will  be  the  0  -  o  direction  (G  =  maximum  in  this  direction)  and  this  is 
the  only  direction  in  which  every  element  of  the  zone  can  migrate  at  the 
same  velocity. 

When  the  atomic  kinetics  are  anisotropic,  the  situation  is  quite 
different  than  the  foregoing.  Anisotropic  p  implies  anisotropic  solid-melt 
surface  energy,  7.  In  fact,  for  any  particular  crystallographic  system, 
a  plane  which  exhibits  a  low  value  of  7  generally  also  exhibits  a  low  value 
of  p.  When  the  anisotropy  of  7  is  sufficient,  stable  microscopic  facet 
formation  will  occur  on  certain  areas  of  the  zone  surface  to  minimize  the 
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free  energy  of  the  system.  Since  the  anisotropy  of  7  is  expected  to 

increase  as  the  temperature  difference  between  the  zone  surface  and  the 
melting  temperatures  of  pure  solvent  increases,  we  can  generally  expect 
a  large  anisotropy  in  p  in  studies  of  alloy  zone  migration.  Thus,  the 
zone  surface  can  be  expected  to  be  composed  of  segments  of  the  "equilibrium 
form,"  the  size  of  the  individual  facets  decreasing  as  G  increases. 

To  decide  whether  the  zone  can  migrate,  with  stable  shape,  in 
a  different  direction  than  the  9  =  o  direction,  it  is  sufficient  to  con-  . 
aider  a  slab  zone  and  a  spherical  zone  moving  in  some  other  0  direction. 
Under  the  assumption  that  the  temperature  distribution  remains  unchanged 
(only  true  for  the  spherical  zone),  VQ  is  given  by  eqs.  12,  13,  and  14 
with  G,  p  and£  replaced  by  G  cos  9,  pQ  and  respectively.  However, 
the  slab  zone  is  easily  seen  to  be  unstable  since  £  is  a  function  of 
and  CQ  is  a  function  of  position  on  the  slab  surface.  Thus,  the  slab  zone 
should  break  up  into  segments .  The  sphere  is  also  unstable  because  the 
isotherms  are  still  normal  to  the  Z  direction  and  must  be  normal  to  the  9 
direction  for  VQ  =  constant  at  all  points  of  the  surface  (V0  =  constant 
is  the  condition  for  zone  shape  stability) . 

Since  the  zone  surface  will  exhibit  facets  of  the  equilibrium 
form  faces  and  since  the  macroscopic  shape  is  unstable  if  the  zone  moves 
in  other  than  the  9  =  o  direction,  the  zone  is  only  stable  if  the  9  =  o 
direction  is  the  outward  direction  -  normal  of  a  face  of  the  equilibrium 
form.  It  must  be  the  most  stable  face  (most  singular) to  assure 
complete  stability.  If  a  cylindrical  or  spherical  zone  is  constrained  to 
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move  in  a  plane,  then  a  similar  stability  condition  applies;  i.e.,  the 
zone  is  stable  if  the  0  direction  is  the  outward  direction  -  normal  of 
one  of  the  lines  of  the  two-dimensional  equilibrium  form  in  this  plane. 

Discussion 

The  experimental,  work  of  Vernick  on  the  migration  of  molten 
Al-Ge  wires  along  a  GO  surface  and  through  a  Ge  block  reveals  a  V  vsi 
relationship  for  this  system  that  is  best  described  by  eq.  13.  This  data, 
which  shows  that  atomic  kinetics  were  important  in  his  experiments,  is 
reproduced  in  Figure  9  along  with  the  best-fit  theoretical  eqs.  13  and  14. 
Eq.  12  was  not  drawn  on- this  graph  since  a  value  of  V  ^  2  x  10 ^  would 
be  needed  for  any  appreciable  fit  and  this  value  is  much  too  large  to  be 
reasonable  for  a  uniform  interface  advance  mechanism.  Reasonable  values 
of .  P  produce  a  V  vs  JL  dependence  only  in  the  range  _JL  ~  10 -  10  ?  cm. 
Likewise  when  one  investigates  the  value  of  A  that  one  might  expect  for 
this  system  using  eq.  4,  one  finds  A  =  500,  a  value  much  larger  than  that 
needed  to  give  a  possible  fit  in  Figure  9*  Reasonable  values  of  A  produce 
a  V  vs  £  dependence  in  the  ranged  ~’10_1  cms. 

The  growth  mechanism  in  these  experiments  seem  definitely  to  be 

a  screw  dislocation  mechanism  with  P  «  88.  The  determination  of  this 

value  of  P  required  the  selection  of  G  =  G*  and  £/m  =  O.75  x  10-^  using 

T  -  66fc,^  1 =  -8°C/At#/13)  C/jy  =  50  At#/13)  kQ  =  2^13)  and  choosing 

-4  2/ 

D  ■  3  x  10  cm  /sec.  This  value  of  D  may  seem  to  be  slightly  larger  than 
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expected..;  however,  it  is  not  inconsistent  with  other  solvent  diffusion 


studies  in  liquid  alloy  systems 


(14) 


We  can  compare  the  measured  value 


of  p  (2/r  )2  =  5  x  10"4  with  the  theoretical  value  from  Hillig  and 

/  ON 

Turnbull's  equation ^  ' 


_  3D  (A  sr 

If-TtyRT.  v 

'  x  m 


(19) 


where  AS  is  the  molar  entropy,  R  is  the  universal  gas  constant,  T.  is 

the  temperature  of  the  system  and  V  is  the  molar  volume.  Eq.  19  gives 

-2 

p  ■sail. 5  x  10  using  the  above  value  of  D.  The  two  results  can  only 

be  reconciled  by  choosing  a  much  smaller  value  of  D  in  eq.  19  than  in 

the  evaluation  of  €  ,  by  assuming  that  G  >  G*  so  that  D  can  be  reduced 

in  £.,  or  by  allowing  that  the  layer  edges  may  be  diffuse.''10^ 

The  effect  of  zone  orientation  on  zone  stability  was  also  in- 
( 1) 

vestigated  by  Wernick.  His  results  are  reproduced  in  Table  I.  In 
this  table,  the  plane  in  which  the  wire  is  constrained  to  move,  the  wire 
direction  in  this  plane  and  the  direction  of  G*  max  (Z-axis)  has  been 
tabulated.  The  direction  of  the  densest  packed  lines  (equilibrium-form 
lines)  are  also  tabulated.  It  may  be  seen  that, except  for  the  (lll):£llo] 
experiment,  the  zone  is  always  unstable  unless  the  9  =  o  direction  was  the 
outward  direction  -  normal  of  the  equilibrium  form.  It  is  interesting  to 
note  that,  for  zone  motion  on  the  (ill)  plane  with  a  jllOj  wire  direction, 
the  zone  should  be  unstable  if  the  9  =  o  direction  is  the  (ll2j  direction 
but  should  be  stable  if  the  9  =  o  direction  is  the  [ll2]  direction.  This 
is  because  the  -;111-  plane  equilibrium  form  is  triangular  rather  than. 


hexagonal  due  to  the  "pseudo-polar’’ -bonding  nature  on  this  plane.  WemickVJ' 
found  that,  in  one  case,  such  a  wire  zone  was  stable  and  this  leads  one 
to  speculate  as  to  whether  the  ©  =  o  direction  was  the  outward  or  the  in¬ 
ward  direction  normal  in  this  particular  experiment. 

Conclusions 

The  migration  rate  of  a  molten  zone  through  a  slab  of  material 
generally  depends  upon  orientation-dependent  atomic  kinetics  of  the  melting 
and  freezing  processes.  The  magnitude  of  the  atomic  kinetic  effect  depends 
upon  the  size  of  the  molten  zone  and  this  size  dependence  of  the  migration 
velocity  may  be  used  as  a  method  for  studying  the  atomic  kinetics  of  alloy 
melting  and  freezing.  The  migration  velocity  is  found  to  be  directly 
proportional  to  the  temperature  gradient  in  the  system.  Stable  zone 
migration  occurs  only  for  migration  in  certain  crystallographic  orienta¬ 
tions  in  a  particular  system.  ' 
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Table  I-The  effect  of  zone  orientation  on  zone  stability  for  Al-Ge  wire  zones 


Plane 

Wire 

Direction 

®max 

Direction 

Most  Dense  Packing 
Directions 

Wire  Stability 

(111) 

[112] 

[iio] 

[Ml],  [oil],  [no] 

Unstable  -  Zone  moves  to  lie  in  [lOl]and  [oil] 
directions 

(111) 

■  [no] 

[H2] 

[Toi],  [oil],  [no] 

Unstable  -  Zone  breaks  up  and  moves  to  lie 
in  <  110>  directions 

(001) 

[010] 

[100] 

[no],  [no] 

Unstable  -  Zone  breaks  up  and  moves  to  lie 
in  <110>  directions 

(001) 

[no] 

[no] 

[no],  [no] 

Stable  -  Zone  travels  as  a  unit  in  [  lio]  direction 

(110) 

[001] 

[no] 

[no] 

Unstable  -  Zone  breaks  up  and  moves  to  lie  in 
<  112>  directions 

(110) 

[no] 

w 

[no] 

Stable  -  Zone  travels  as  a  unit  in  the  [ooi] 
direction 
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FIGURE  CAPTIONS 


Figure  1 

Figure  2 
Figure  3 

Figure  4 

Figure  5 

Figure  6 

Figure  7 

Figure  8 

Figure  9 


„  Portion  of  phase  diagram  plus  the  solid  and  liquid  charge 
used  in  temperature  gradient  zone  melting. 

.  Spatial  orientation  of  slab,,  cylindrical  or  spherical  zone. 

.  Portion  of  phase  diagram  illustrating  both  the  equilibrium 
and  actual  temperature  and  solvent  concentration  at  the  zone 
surfaces  for  Q  -  0  and  ©  ~7T- 

.  Thermal  parameter  plot  for  different  zone  shapes  allowing  the 
determination  of  G  from  a  knowledge  of  the  thermal  parameters 
of  the  system,  and  the  gross  experimental  conditions . 

.  Plot  of  the  relative  zone-edge  migration  rate  -V/v  o 
function  of  the  dimensionless  concentration  parameter 
(C^-CqJ/Ctp  for  several  values  of  the  solvent  partition 

coefficient,  k  . 

3  o 

.  Plot  of  zone  migration  velocity*  V<^=*  as  a  function  of  zone 

size*  JL  *  for  the  diffuse  interface  mechanisms  (T’s  l/u^»-  l/p 
and  the  screw  dislocation  mechanism  »  l//py,+  1  /vPq*  h 

.  Plot  of  the  zone  migration  velocity*  V*  as  a  function  of  zone 
size*  &-  *  for  the  two-dimensional  nueleation  mechanism  assuming 
Ao/ln  (d0  / V-jy.  )<.<kjf~j In  (n  -jj/ V -j?-)  « 

.  Solvent  distribution  in  the  zone*  C*  as  a  function  of  distance* 
Z*  for  the  conditions  of  (i)  infinite  atomic  kinetics  with  no 
mixing  and  (ii)  finite  atomic  kinetics  with  no  mixing*  partial 
mixing  and  complete  mixing. 

(■*) 

.  Fitting  of  Wernick’s  dataw/  on  the  migration  of  Al-Ge  wires 
through  Ge  with  eq.  13  and  eq.  14  choosing  best-fit  parameters. 


(a) 

Slab 


(b) 

Cylinder 


Fi«.  2 


(c) 

Sphere 
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